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Abstract 

We discuss a scenario that gravitinos produced non-thermally by an inflaton decay constitute 
dark matter in the present universe. We find that this scenario is realized for wide ranges of the 
inflaton mass and the vacuum expectation value. What is intriguing about this scenario is that 
the gravitino dark matter can have a relatively large free streaming length at matter-radiation 
equality, which can be probed by future observation on QSO-galaxy strong lens system. 
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In spite of accumulating observational data supporting the presence of the dark matter 
(DM) in our universe we have not yet identified what DM is made of. Among many 
candidates proposed thus far, the gravitino, a supersymmetric parter of the graviton, is 
particularly interesting, and it has been thoroughly investigated in connection with leptoge- 
nesis [2j and the collider signatures [3|. 

The gravitinos are copiously produced by particle scattering in thermal plasma, once the 
decay of the inflaton reheats the universe. If the gravitino is the lightest supersymmet- 
ric particle (LSP), it is stable and can be a good candidate for DM The 
gravitino abundance is directly related to the reheating temperature, Tr. In particular, for 
the gravitino mass m 3 / 2 larger than O(10)GeV, the required reheating temperature for the 
gravitinos to be DM is so high, T R > 10 9 GeV P, Q, that the thermal leptogenesis scenario 
may work [ll]. 

However, the detailed studies on the big bang nucleosynthesis (BBN) revealed that the 
abundance and the ^lifetime of the next-to-lightest supersymmetric particle (NLSP) are 
tightly constrained [12']. This drove the above attractive scenario into a corner, since the 
lifetime of NLSP tends to be longer than the BBN bound especially for 777,3/2 larger than 
O(10)GeV, signaling the need for some changes. Several solutions has been proposed: e.g ., 
a late-time entropy production HQ , [ll]] and a theory with R-parity violation [ly, u\ - 



Another is to abandon thermal leptogenesis and consider a non-thermal leptogenesis sce- 
nario instead, which requires a lower reheating temperature, Tr > 10 6 GeV. Then 
the gravitino can account for the observed relic density even for 777.3/2 lighter than O(10) GeV, 
making it easier to evade the constraints from the NLSP decay. One drawback of this ap- 
proach however is that one needs to introduce ad hoc couplings of the inflaton with the 
right-handed neutrinos. 



Furthermore, it has been recent 
by an inflaton decay 
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y pointed out that the gravitinos are generically produced 



26(| . Since such non-thermal gravitino production 



generically occurs for the most inflation models, it is worth studying how it affects the 
conventional picture on the gravitino DM scenario. In this letter, we pursue a possibility 
that the gravitinos produced by an inflaton decay constitute a dominant component of DM. 
Generically, the required reheating temperature becomes lower than without the non-thermal 
production. This makes it difficult to integrate the thermal leptogenesis scenario into this 
framework. As we will see later, however, the right-handed (s)neutrinos are generated by 
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the inflaton decay, and their subsequent decay may generate a right amount of the baryon 
asymmetry via leptogenesis for certain values of the inflaton parameters [2^]. What is 
particularly appealing about this scenario is that both the gravitino DM and the non- 
thermal leptogenesis can be realized without introducing any couplings ad hoc by hand, 
if the inflaton parameters satisfy certain conditions. In addition, the produced gravitinos 
can have a large velocity at matter-radiation equality, which affects the growth of density 
fluctuations of DM a . Future observations on e.g. QSO-galaxy strong lens system {32] may 
be able to support or refute this scenario. 

Let us first briefly review the recent development on the gravitino production from the 
inflaton decav There are three gravitino production processes; (a) the gravitino pair pro- 



duction 
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23|; (b) spontaneous decay at tree level |24|; (c) anomaly- induced decay 
at one-loop level 25]. For the processes listed above, the gravitino production rate can be 
expressed as 

r x (WY m * m 

Fs/2 " 32^ J M*' { ) 

where is the inflaton mass, ((f)) a vacuum expectation value (VEV) of the inflaton, and 
Mp = 2.4 x 10 18 GeV the reduced Planck mass. Here it should be noted that (<fi) is evaluated 
at the potential minimum after inflation. The precise value of the numerical coefficient x 
depends on the production processes, possible non-renormalizable couplings in the Kahler 
potential, and the detailed structure of the supersymmetry (SUSY) breaking sector [26;]. To 
be concrete, let us assume the minimal Kahler potential and the dynamical SUSY breaking 



(DSB) 33] with a dynamical scale A. In the DSB scenario, the SUSY breaking field z 



can acquire a large mass m z , which is assumed to be roughly equal to the dynamical scale 



A ~ ym3/2Mp in the following. Such a simplification does not essentially change our 
arguments. For a low-inflation model with < A, the process (a) becomes effective, and 
x — 1. On the other hand, for the inflaton mass larger than A, the processes (b) and (c) 
become effective instead. The inflaton decays into the hidden quarks in the SUSY breaking 
sector via Yukawa couplings (process (b)), or into the hidden gauge sector via anomalies 
(process (c)). Since the hidden quarks and gauge bosons (and gauginos) are energetic when 
they are produced, they are expected to form jets and produce hidden hadrons through the 



Such a DM candidate with a large velocity and its astro-physical implication was first discussed in 
Refs . 



28l I29L |30| , and intensively studied in connection with the so-called super WIMP mechanism [3l| 
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strong gauge interactions. The gravitinos are likely generated by the decays of the hidden 
hadrons as well as in the cascade decay processes in jets. We denote the averaged number 
of the gravitinos produced per each jet as i\T 3 / 2 . Then x is given by [26] 



b 



■'■ ^(\N y \YZ\+N g al{T^-T^Y), (2) 



in 2 



where Yh and ah are the Yukawa coupling and a fine structure constant of the hidden gauge 
group, respectively, N y denotes a number of the final states for the process (b), N g is a 
number of the generators of the gauge group, and and TjP are the Dynkin indices of 
the adjoint representation and the matter fields in the representation r. Although x depends 
on the structure of the SUSY breaking sector, its typical magnitude is 0(1CT 3 — 1CT 2 ) for 
777.0 > A c . To be concrete we take 

«=( 1 fOT m * <A , (3) 

I 10 3 or 10~ 2 for m 4> > k 

in the following. 

Using the gravitino production rate given above, we can estimate the abundance of the 
gravitinos non-thermally produced by an inflaton decay: 



Y 



(NT) _ T 3 /2 

" v2007 ^10 15 GevJ Uo 12 GevMl0 6 GeVy ' [ ) 

where g* counts the relativistic degrees of freedom, and denotes the total decay rate of 
the inflaton that is related to the reheating temperature as 

Equivalently, the gravitino density parameter is 

°SP* - «»»* &) (lo^v) 2 (lo^v) 2 (JkvT ■ W 



b If the Kahler potential takes a form of the sequestered type, the spontaneous decay through Yukawa 
couplings is suppressed [2J, |25 1 . 



Roughly we expect N 3/2 = 0(1 - 10 2 ), N g = 0(1), a h = 0(0. 1), and T^ h) - T^ h) = 0(1), while Y h 
strongly depends on the SUSY breaking models. Note also that the gravitino can be produced through 
the Yukawa interaction in the messenger sector, if the inflaton mass is larger than the messenger scale. 
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where h is the present Hubble parameter in units of lOOkm/s/Mpc. Note that the gravitino 
abundance is inversely proportional to the reheating temperature. Due to this feature, the 
non-thermal gravitino production tends to require a relatively low Tr to realize the gravitino 



DM scenario. Indeed, by solving VL^J^h 2 = 0.11 



"3/2 



with respect to T R , we obtain 



So, if Tr is given by the above value, the non-thermally produced gravitino has a right 
abundance to become a dominant component of DM. Generically, one has to introduce a 
coupling of the inflaton to the standard-model sector with an appropriate strength, in order 
to realize Tr given by Eq. (J7J) . However, there is a natural way to induce the reheating, and 
we will discuss this possibility later. 

For the non-thermally produced gravitinos to account for DM, several conditions must 
be met. First, the gravitino production by thermal scatterings should give only negligible 
contribution to the D1V abundance. The abundance of the gravitinos produced by thermal 



35] 



scatterings is given by [y, 

n^h 2 ~ o 14 ( m ' 93 V f^H-Y 1 ( Tr ^ r 8 ) 

3 / 2 " U v300GeW llGeW llO 8 GeW ' [ ) 

where m~ gz is the gluino running mass evaluated at T = Tr. Requiring Q^h 2 to be less 
than the observed DM abundance, Qdm^ 2 — 0.11, T R is bounded above: 

This constraint is valid for m 3 / 2 > 100 keV, which is satisfied for the parameter space 
concerned as shown later. 

Another constraint comes from the recent discovery that, once the inflaton acquires a 
non-vanishing VEV, the inflaton decays into the visible sector through the top Yukawa 



coupling 241- Due to the presence of this decay process, Tr cannot be arbitrarily low. 
Indeed, it is bounded below as 

^^^^(A)^)' do, 

where Y t is the top Yukawa coupling. The inequality is saturated if the inflaton has no direct 
couplings with any other fields in superpotential d . 



Note that we assume the minimal Kahler potential in the Einstein frame. 
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The last constraint arises from the fact that the non-thermally produced gravitinos can 
have a large velocity at matter-radiation equality, in contrast to the gravitinos produced by 
thermal scatterings. This not only limits the parameter space, but also provides a possibility 
that the scenario may be probed by future observation on QSO-galaxy strong lens system. 
Let us estimate the comoving free streaming length of the gravitino at matter-radiation 
equality, assuming that it has an initial energy, when produced. For < A, we 

have e = 1 since a pair of the gravitinos is directly produced by the inflaton decay. On the 
other hand, for > A, multiple gravitinos are indirectly generated by the inflaton decay, 
and so, its energy tends to be smaller than we expect e ^ 0{N^ 2 ) = 0(1CT 3 — 0.1). 

To be concrete we will take e = 10~ 3 or 10 -2 for > A. The comoving free streaming 
length Xps at matter-radiation equality is defined by 

Jt D a(t) 

where a(t) is the scale factor, and to and t e q(~ 2 x 10 12 sec) denote the time at the inflaton 
decay and at matter-radiation equality, respectively. t> 3 / 2 is the velocity of the gravitino, 
given by 

IPS/2 





fe)) 


2 





where we have approximated 3> m 3 / 2 , and a D is the scale factor at the inflaton decay. 
Integrating ffTTj) yields 

1 



-.X- 1 sinh" 1 X, 



H 0y /l + z eq 

-*-<-»(l) 4 (S)li^v)(iolv) _1 W) 



with 



_ 2m 3/2 a eq 

A = , 

erritf, a D 

~ 8 x 10 2 6- 1 f^-V 4 ( m + Y 1 ( Tr ) (14) 

" \2007 \lGeVj ll0 12 GeW Uo 5 GeW ' K > 

where Ho is the Hubble parameter at present, and z eq and a eq are the red-shift and the scale 
factor at the matter-radiation equality. In the second equation of ([TBI , we have assumed 
X ^> 1 and used Hq 1 ~ 4 x 10 3 Mpc and z eq ~ 3000. In Eq. (|T4|) . we have used ao/a eq = 
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(r0-t eq ) 1 / 2 : . The constraint from Ly-a clouds, Xps ^ 1 Mpc, implies X > 450. We therefore 
obtain a constraint on Tr as 

rs > 5xl0lGeV£ (^)-'(^)-'(_^_). (15) 

The meaning of this constraint is clear: the reheating must occur so early that the velocity 
of the produced gravitino becomes small enough due to redshift by the matter-radiation 
equality. 

Thus, if the reheating temperature Tr is given by ([7]) and satisfies the above constraints 



©, ([TOD, and (USD in addition to the BBN constraint T R > lOMeV [33, |3J, |39|, the non- 
thermally produced gravitinos account for DM. As mentioned above, one may have to add 
appropriate couplings of the inflaton to light degrees of freedom, in order to realize Tr 
given by (j7]). However there is one interesting possibility that the reheating is induced by 
the decay through the top Yukawa coupling. Then the inequality (flO]) becomes saturated. 
This is the case if there are no direct couplings of the inflaton with any other fields in the 
superpotential. The presence of the decay process through the top Yukawa coupling not 
only constrains the reheating temperature, but also provides an intriguing way to induce 
the reheating. For the moment let us pursue this possibility. From (J7J) and ( ITUl) . we obtain 



W ^ ( m * )* ~ o 01 IKIz- 1 (-^W-^ 
10 15 GevHl0 12 GeW " v2007 llGeV 



(16) 



Thus, if the inflaton parameters, and (0), satisfy the above relation (|16p . the non- 
thermally produced gravitino has a just right abundance to be DM. Interestingly, the free 
streaming length becomes independent of the inflaton parameters and the gravitino mass in 
this case. Indeed, Xps is approximately given by 

^--M^Wfe)- (17) 

The Ly-a constraint requires ex < 10~ 4 , which is naturally satisfied for a high-scale inflation 
model with m§ > A. It is intriguing that the gravitino DM scenario points to a high- 
scale inflation model with > A and predict a relatively large free streaming length, as 
long as the reheating is induced by the top Yukawa coupling. For x = 10~ 3 ~ 10~ 2 and 
e = 10 -3 ~ lO^ 1 , the comoving free streaming length takes a value from 10 kpc up to 1 Mpc 
(limited by the Ly-a constraint). 
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Now let us consider the inflaton decay into the right-handed (s)neutrinos thorough large 
Majorana mass terms: 



(18) 



where i = 1,2,3 is the family index. We consider the inflaton decay into the lightest right- 
handed (s)neutrino Ni for simplicity, assuming that the decay into the heavier ones, N2 and 
iV 3 , are kinematically forbidden. We drop the family index in the following. The partial 
decay rate of the inflaton into the right-handed (s) neutrinos is [cf. 241 ]] 

2 



1 



A' 



16vr \Mp 



M 2 P 



\ 



1 



4M 2 



m 



2 ' 



(19) 



where we have taken account of both the decay into the right-handed neutrinos and that into 
the right-handed sneutrinos. Note that one does not have to introduce any direct couplings 
of the inflaton with the right-handed neutrinos to induce the decay. The decay proceeds as 
long as the inflaton acquires a nonzero VEV. 

The lepton asymmetry can be produced by the decay of the right-handed (s) neutrinos, if 
CP is violated in the neutrino Yukawa matrix |2j. The resultant lepton asymmetry is given 
by 

(20) 



n L 3 T R 

S 2 771$ 



where Bn = IV /T^ denotes the branchingratio of the inflaton decay into the (s)neutrinos. 



The asymmetry parameter e\ is given by 



4Q| 



(21) 



where m V3 is the heaviest neutrino mass and <5 e ff < 1 represents the effective CP- violating 



phase. The baryon asymmetry is obtained via the sphaleron effect: 



4l| 



n B 
s 



23 s 



(22) 



Using the above relations, we obtain the right amount of baryon asymmetry, 

3 / ija \ 2 , m x -1 



* lxlO- 9 (^-Y~ 2 
s V200/ 

~ 5xl0- n f-^f* 
V200/ 



M 



10 13 GeV 



(0) 



T, 



R 



10 16 GeV/ V10 6 GeV 



M 



(0) 



10 13 GeV/ \10 16 GeV/ V10 14 GeV 



m<t> 



771,, 



0.05eV 



$ e g, (23) 
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where we have assumed that the inequality (fTD|) is saturated in the second equality. Note 
that M cannot exceed Therefore, the baryon asymmetry is proportional to positive 

powers of m^, if M is set to be a value that maximizes the asymmetry. 

In Figs. [T] and [21 we show the parameter space where the reheating temperature ([7]) 
satisfies the above constraints 0, (HOI) , and (1151) . in addition to the BBN constraint Tr > 
10 MeV. In the shaded (blue) regions, the baryon asymmetry can be explained by the non- 
thermal leptogenesis scenario discussed above, if an appropriate value of M(< 10 15 GeV) is 
chosen (we set m Vz = 0.05 eV and 5 e fr = 1). We have chosen several values of the gravitino 
mass: 1713/2 = 100 MeV, 1 GeV, and 10 GeV e . For smaller 1713/2, one needs to generate more 
gravitinos, due to which the allowed region shifts upward. At the same time, the constraints 
from and f[T5"|) become severer for smaller m 3 / 2 , reducing the allowed space. Thus, if 
x = 10~ 3 (10~ 2 ) for m<f> > A, the gravitino mass should be larger than 1 MeV (100 keV) 
for the non-thermally produced gravitinos to account for DM, since otherwise there is no 
allowed region for (0) < M P . If x becomes larger for > A, the allowed region shifts 
downward, and a smaller value of the gravitinos mass becomes allowed. On the other hand, 



if x becomes smaller due to e.g. conformal sequestering 26j, we have more parameter space 
for the non-thermal leptogenesis to work successfully. 

The dotted (red) lines correspond to the special case that the reheating is solely induced 
by the decay through the top Yukawa coupling and the non-thermally produced gravitinos 
become DM. Therefore the inflaton parameters on the dotted (red) lines are particularly 
interesting in a sense that one does not have to introduce any couplings ad hoc by hand; the 
decay spontaneous proceeds through the top Yukawa coupling, and the gravitino has just a 
right abundance to become DM. Note that the free streaming length is constant along the 
dotted (red) lines and independent of 717,3/2, and (</>), as mentioned before. 

In Figs. [1] and [21 we also show the contours of the free streaming length Xps — 1 kpc, 
10 k pc, 100 kpc, and 1 Mpc. The future submillilensing observations can cover Xps ~ 2 
kpc [32]. In particular, since the interesting case that the reheating occurs through the 
top Yukawa coupling (dotted red lines) predicts the gravitino DM with a relatively large 
free streaming length (> 10 kpc), it can be probed by future observations. Such a large free 



Note that we have not taken account of the constraints on the NSLP decay. Therefore, the figure in the 
case of m 3 / 2 = 10 GeV is valid only if the cosmological problems associated with the NLSP are somehow 
avoided by e.g. introducing R-parity violating operators with an appropriate magnitude [l^. [17I]. 
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streaming length may also solve the missing satellite problem [45| and the cusp problem [46 ]. 

From the figures, one can see that relatively broad ranges of the inflaton mass and VEV 
are allowed. In particular, when combined with the non-thermal leptogenesis scenario, we 



are led to a high-scale inflation model with > A. However, studying 



spaces of the representative high-scale inflation models (such as the hybrid 42| and smooth 



the parameter 



hybrid [43(, and chaotic 44] f inflation models) in detail, one finds that only small part of 
the parameter space actually overlaps with the region where the non-thermal leptogenesis 
works, especially if x takes a value on the high side ~ 10~ 2 ; those inflation models tend to 
predict lighter m§ and larger (0) compared to those favored by the non-thermal leptogenesis. 
See Fig. [3J Such a tension may be ameliorated if one assumes some mechanism (e.g. the 
conformal sequestering) to suppress x to a smaller value. 

In summary, we have considered a scenario that the non-thermally produced gravitinos 
from the inflaton decay become a dominant component of DM. Interestingly, if the reheating 
is induced solely by the decay through the top Yukawa coupling, a high-scale inflation 
model is required for the non-thermally produced gravitinos to account for DM, and the free 
streaming length Xps is predicted to be in the range between O(10) kpc and O(0.1)Mpc, 
independently of the inflaton parameters and the gravitino mass. Such large free streaming 
length may affect the growth of the density fluctuations in DM. The suppression of the 
density contrast below the free streaming scale results in the absence of the sub-halos. This 
feature may be supported or refuted by future observations on the QSO-galaxy strong lens 



system [32]. 
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m 3 /2 = 10 GeV 




[GeV] 

FIG. 1: In the regions surrounded by the solid (green) lines, the non-thermally produced gravitinos 
can account for the observed DM density, if Tr is given by ([7]). We have imposed the constraints 
from (i) thermal production of the gravitino (see ([9])); (ii) decay through the top Yukawa coupling 
(see (fTUj) ): (hi) Ly-a clouds (see (fT5j) ). On the dotted (red) line, the reheating is induced solely by 
the decay via the top Yukawa coupling and the non-thermally produced gravitino explains DM. 
The thin solid (orange) lines are the contours of the free streaming length Xps = 1 kpc, 10 kpc, 100 
kpc, and 1 Mpc, from top to bottom. In the shaded (blue) regions the present baryon asymmetry 
can be explained by the non-thermal leptogenesis. We set x = 10~ 3 and e = 10~ 2 for > A, 
respectively. 
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FIG. 2: Same as Fig. [T] except for x = 10 2 and e = 10 3 for > A. 
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FIG. 3: We show the representative high-scale inflation models; hybrid 
(blue) line), smooth hybrid 



42j ] (thick long dashed 



441 ] (long dashed 



431 ] (thick short dashed (purple) line), and chaotic 
dotted (red)) inflation models, superposed on the panels of m 3 / 2 = lGeV shown in Fig. Q] (left) 
and Fig. [2] (right). 
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